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Abstract. Sunspot fine structure has been modeled in the past by a combination of idealized 
magneto-convection simulations and simplified models that prescribe the magnetic field and 
flow structure to a large degree. Advancement in numerical methods and computing power has 
enabled recently 3D radiative MHD simulations of entire sunspots with sufficient resolution to 
^ , address details of umbral dots and penumbral filaments. After a brief review of recent develop- 
ments we focus on the magneto-convective processes responsible for the complicated magnetic 
structure of the penumbra and the mechanisms leading to the driving of strong horizontal 
outflows in the penumbra (Evershed effect). The bulk of energy and mass is transported on 
scales smaller than the radial extent of the penumbra. Strong horizontal outflows in the sunspot 
penumbra result from a redistribution of kinetic energy preferring flows along the filaments. This 
redistribution is facilitated primarily through the Lorentz force, while horizontal pressure gra- 
dients play only a minor role. The Evershed flow is strongly magnetized: While we see a strong 
reduction of the vertical field, the horizontal field component is enhanced within filaments. 
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1. Introduction 

Sun and starspots play a central role for our understanding of solar and stellar mag- 
netism. Since direct observations of magnetic field in stellar convection zones are very 
limited (e.g. helioseismic inversions), sun and starspots provide a window to understand 

>> \ ' the magnetism of stellar interiors, provided we understand in detail how starspots form, 
evolve and decay. Sunspots are a multi scale problem with respect to spatial and tempo- 

y—i ■ ral scales. While their typical size is > 20 Mm, fine structure is observed at the smallest 
scales currently obs ervable of about 200 k m (e.g. high resolution ground based observa- 
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sr<! e.g. Ilchimoto et al. ( 20071 )). Properly resolving the scales currently observed requires 



tions w ith the SST, Scharm er et al.l ([20021 ). or space based observations with HINODE, 



grid resolutions of 20 km or less. Details of the penumbra evolve over time scales of hours, 
while the life time of sunspots and the evolution of the adjacent moat happens on time 
scales of days to weeks. On the other hand typical numerical time steps are of the order 
of 0.1 sec (assuming already that the very fast Alfven velo cities found above a s unspot 
umbra were removed from the system, see e.g. appendix of Rempel et al. ( 2009al )). As a 



consequence a well resolved realistic numerical simulation of an entire sunspot requires 
billions of grid points and hundred thousands of time steps. Since the essential physics 
that need to be considered (MHD, 3D radiative transfer, realistic equation of state) are 
well known and have been included in MHD codes for more than 2 decades, the primary 
challenge for addressing sunspot structure is linked to robust efficient numerical schemes 
and availability of powerful computing resources. The latter became finally available on 



1 



2 



Matthias Rempcl 



the scale needed, which allowed for substantial progress within a time frame of only a 
few years. In the following sections we will very briefly summarize recent progress and 
a few aspects of sunspot structure for which numerical simulations provided substantial 
insight. 



2. Recent developments 

We limit the following discussion entirely to MHD simulations that include 3D ra- 
diative transfer and a realistic equati on of state. Here, p r ogres s in applications related 
to sunspots started with the work bv lSchussler k, Vogler (2006) who presented a MHD 
simulation of a sunspot umbra showing the development of magneto-convection in form 
of umbral dots. The simulation revealed that almost field free upflow plumes can form 
withi n an initially monolithic umbra and transport energy through overturning convec- 
tion. iHeinemann et al" ( 2007 ) focused on a narrow slab through the center of sunspot. 
Their setup allowed to model the transition from umbra toward penumbra and granu- 




Figure 1. Numerical sunspot model. The domain size is 49 x 49 x 6.1 Mm, the grid resolution 
16 x 16 x 12 km (3072 x 3072 x 512 grid points) . The simulation was performed with grey radiative 
transfer, displayed is the bolometric intensity for a vertical ray in the range from 0.25 — 1.5/0. 
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lation, while keeping the computational expense at a moderate level. Their simulation 
showed the formation of short filaments with dark lanes at the umbra/granulation in- 
terface, however, umbral dots were not present. The filaments were propagating inward 
duri ng the formation phas e and showed outflow along their axis. Based on this simula- 
tion Scharmer et al.l (|2008f) interpreted the Evershed flow as convective flow component 



along filaments. In a, very similar setup with an initially larger segment of a sunspot 
Rempel et al. (2009a) were able to produce umbral dots and filaments of about 3 Mm 



length, showing the smooth transition from central to peripheral umbral dots and fila- 
ments in the inner penumbra. The latter show the presence of a bright head propagating 
inward and detaching from the filament during the formation phase and up to 3 Mm long 
dark lanes with evidence of twisting motions. Along filaments, outflows of a few km/s were 



presen t, and on a larger scale a moat flow, previously also reported in lHeinemann et al 



( 20071) . Overall these simulations clearly stress the common magneto-convective origin 
of umbral dots and penu mbral filaments, but a realistic outer penumbra with strong 
outflows was not present. Kitiashvili et al.l (2009) conducted a study of how the mag- 



netic field strength and inclination angle influence the formation of horizontal flows in 
magneto-convection. Fast flows with mean amplitudes of I — 2 km/sec were found for 
a 1.5 kG field with inclination angle of 85 degrees. The fir st comprehensive simulation 
of entire sunspots was presented bv iRempel et al. ( 2009b ). Using a setup including a 



pair of opposite polarity sunspots in a 98 x 49 Mm wide and 6.1 Mm deep domain, this 
simulation showed the formation of extended (up to 10 Mm wide) penumbrae with hor- 
izontal outflows of up to 6 km/sec mean- and 14 km/sec peak flow speeds. At the same 
time this simulation contained umbral dots and a substantial moat region, connecting 
all the aspects of sunspot fine structure in one comprehensive simulation run. While the 
penumbra con t ained r adially elongated convection cells, the intensity image (see Fig. 1 in 
Rempel et ail (|2009bl )) did not yet show the typical radial structure of narrow filaments 



known from observatio ns. The currently bes t resolved sunspot simulation is presented 
in Fig. [I] Compared to IRempel et al. | (|2009bl) the horizontal resolution is doubled (from 



32 to 16 km), while the vertical resolution is increased from 16 to 12 km, leading to a 
total grid size of 3072 x 3072 x 512 for a domain of the size 49 x 49 x 6.1 Mm. Here 
we focused again on a single spot, but artificially enhanced the field inclination through 
the top boundary condition to generate conditions tha t are co mparable to the region in 
between the opposite polarity spots of IRempel et al. ( 2009b ). While the photosphcric 



appearance of the sunspot is substantially improved compared to previous simulations, 
we did not find fundamental differences in the underlying magneto-convection process. 

In the following section we highlight a few central aspects we learned from the simu- 
lations summarized above. 



3. Magnetic fine structure of penumbra 

Figure [2] summarizes the magnetic fine structure of umbra and penumbra at the 
tross = 1 level for the simulation presented in Fig. [TJ A common element in umbra 
and penumbra is a strong reduction of the vertical magnetic field component, which 
is associated with umbral dots, peripheral umbral dots and penumbral filaments. The 
horizontal magnetic field component shows a different behavior. A reduction of the hor- 
izontal field strength is only present in the umbra, whilst in the penumbra filaments 
(flow channels) have an enhanced horizontal magnetic field strength throughout, which 
is evident from the correlations shown in the bottom right panel of Fig. [5] The com- 
bination of reduced vertical and enhanced horizontal field leads to the observationally 
inferred interlocking comb structure with strong variation of the field inclination angle 
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and st rong horizontal flows in the component with nearly horizontal field (see e.g. lSolanki 
( 2003^ The combined effect of the weakening of B z and strengthening of Br within 
flow channels results in an overall reduction of \B\ in the inner, and increase of \B\ in 



the outer, penumbra. Recently llchimoto et al.l (|2OO70 found this trend in the observed 
Vr — \B\ correlation, even though positive values of the correlation were not observed. 
An Evershed flow in regions with enhanced field strength in the outer penumbra was 
proposed bv lTritschler et al.l (|2007h based on zero-crossings of the NCP. 

Overall the numerical simulation presented here supports a strongly magnetized penum- 
bra and Evershed flow. While the underlyin g structure of the magneto -con vection is well 

captured by the "gappy" penumbra model of lScharmer fc Spruit ( 20061 ) and Spruit fc Scharmer 
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Figure 2. Magnetic fine structure of sunspot at r = 1 level. The top panels show vertical 
(±3 kG) and radial field strength (0 ... 3 kG). The bottom left panel shows the inclination, with 
black indicating vertical outward, grey horizontal, and white vertical inward directed field. We 
set the inclination to for regions with less than 500 G field strength (dark color outside the 
penumbra). The bottom right shows correlations between horizontal flow velocity and magnetic 
field strength fluctuations as a function of the spot radius. The vertical dotted lines correspond 
to the dark concentric circles shown in the other three panels. 
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( 20061) . we do not see any support for the claim that these gaps are close to field free. 



While the latter is a good approximation for the vertical field component alone, h o rizon - 
tal field is actually enhanced compared to the background field. Recently iRempel ( 2010l ) 



showed that the horizontal field originates from the strong subsurface shear profile of 
the Evershed flow, which leads to an induction term B z 8 z vr with a strength of about 
5 — fO G/sec a few 100 km beneath r = 1. Relative contributions from numerical magnetic 
diffusivity are only of a few percent. A recent convergence study (Rempel, in preparation) 
from 96 x 32 up to 16 x 12 km resolution (Fig. [T] presents the highest resolution case) 
strongly supports the robustness of this result. 



4. Convective energy transport 

All numerical simulations performed to date (most of them are summarized in Sect. [2J 
point toward a common magneto-convective origin of sunspot fine structure. Comparing 
the azimuthally averaged bolometric intensity and vertical RMS velocity at the tr oss = 1 
level we find a relationship of the form / oc y/v^- MS (T — 1), which holds very well from 
umbral dots to penumbral filaments and even granulation. From this we can conclude 
that the vertical RMS velocity in the inner penumbra with / w 0.7/0 should be about 
half of the convective RMS velocity seen in granulation, i.e. it should be about 1 km/sec 
instead of 2 km/sec. This is, at least in a s tatistical sense, compatible with the findings 
reported in iFranz fc Schli (2009), who found about 500 m/s vertical RMS in 



the penumbra vs about 1 km/sec RMS in the quiet sun. The shortfall by a factor of 2 
for both penumbra and quiet sun is due to the limited resolution of observations and the 
fact that the typically used spectral lines form at levels higher than tr OS s = 1- 

We find throughout the penumbra an upflow filling factor in the 40 — 60% range and 
most of the mass flux is turning over on scales significantly shorter than the radial extent 
of the penumbra. Only about 15% of the total overturning mass flux (similarly also 
energy flux) is found in the azimuthal mean component (corresponding to an upflow of 
about 200 m/s in the inner and downflow of up to 500 m/s in the outer penumbra). The 
mass flux of the penumbra is balanced within the bounds of the penumbra. 

The almost equal presence of up and downflows everywhere in the penumbra is cur- 
rently not inferred from most observations (which see a preference of upflows in the inner 
and downflows in the outer penumbra) and might require even higher resolution. With- 
out overturning motions on scales much shorter than the radial extent of the penumbra 
it would not be possible to maintain the observed brightness of about 0.7/0. 



5. Origin of Evershed flow 

Recently iRempell (l2010h analyzed in detail the processes underlying the driving of 



the Evershed flow. To this end the contributions from acceleration, pressure, buoyancy 
and Lorentz forces in the kinetic energy equation were compared between the plage 
region (surrounding the sunspots) and the sunspot penumbra. In the plage region most 
of the pressure/buoyancy driving takes place in downflows and is in balance with vertical 
acceleration forces. Horizontal flows are driven by horizontal pressure gradients. In the 
penumbra the pressure/buoyancy driving is shifted into upflow regions and is there in 
balance with vertical Lorentz forces, while acceleration forces are unimportant. In the 
horizontal direction acceleration forces are in balance with the horizontal Lorentz force 
component, while horizontal pressure gradients do not play a major role throughout most 
of the penumbra. The Lorentz force facilitates the energy exchange between motions in 
the vertical and horizontal direction, while the net work done by Lorentz forces remains 
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negative (sink of kinetic energy). Some aspects of this picture (pressure driving in upflows, 
deflection and guiding of motions by magnetic field) have be e n cap t ured to some extent in 
simplifi ed fl ux tube models such a s Montesinos fc Thomas! ( 19971 ). Schlichenmaier et al. 
(1998b) and lSchlichenmaier et al.l ( 1998a ). However, the presence of vigorous convection 
leads to a situation in which the outflowing mass is continuously replaced through upflows 
extending along filaments - a situation that is inherently opposite to the concept of a flux 



tube. Overall the Evershed flow is best characterized as convective flow (jScharmer et al 
although, compared to field free convection, notable differences exist with respect 
to the underlying driving forces. 



6. Conclusions 

Numerical simulations have been advanced to the point where they provide a con- 
sistent and unifying picture of the magneto-convective processes underlying the energy 
transport, magnetic fine structure and origin of large scale flows in sunspots. The bulk 
of energy and mass is transported by overturning convection with scales substantially 
shorter than the radial extent of the penumbra. The filamentation of the penumbra 
and the driving of large scale outflows is strongly linked to the presence of anisotropic 
convective flows. We find the above picture being converged with respect to numerical 
resolution. We explored the range from 96 x 32 km to 16 x 12 km and found that most 
aspects are already well described with 48 x 24 km resolution. For the currently accessible 
resolution range, flows within filaments are mostly laminar; whether a possible transition 
to turbulent flows at higher resolution could change results remains an open question (due 
to the rather strong magnetic field turbulence might remain suppressed even at higher 
resolutions than currently affordable). Initial state and boundary conditions at top and 
bottom have a strong influence on the global magnetic structure and stability (i.e. they 
can make the difference between having or not having a penumbra), however the details 
of sunspot fine structure as discussed here are influenced to a much lesser degree. While 
the overall radial extent of the penumbra is subject to boundary conditions, the details 
of filamentation, energy transport, and driving mechanism of the Evershed flow are not. 
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